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We present a finite-temperature theory of the anisotropic spin-spin correlations in magnetic 
metallic monolayers, deposited on a suitable substrate. The 'spins' are the local moments set up by 
the itinerant electrons, and the key concept is the relativistic disordered local moment state, which 
represents the paramagnetic state of a set of local moments. The spin-spin correlations between 
these local moments are then extracted using the linear response formalism. The anisotropy is 
included in a fully relativistic treatment, based on the Dirac equation, and has a qualitative impact 
on noncollinear magnetic states, by lifting their chiral degeneracy. The theory is applied to Mn 
monolayers on the hexagonal (111) surfaces of Pd, Pt, Ag and Au. The presence of competing 
exchange interactions is highlighted by choosing different substrates, which favour either the row- 
wise antiferromagnetic state or the chiral triangular Neel state. We correlate the electronic structure 
with the magnetic properties, by comparing filled with partially filled substrate d-bands, and low 
vs high atomic number. The disagreement between theory and experiment for Mni/Ag(lll) is 
addressed, and the nature of the magnetic domains found experimentally is suggested to be chiral. 

PACS numbers: 75.70.Ak, 71. 15. Mb, 75.30. Gw, 75.50.Ee 



I. INTRODUCTION 

A thorough theoretical description of the magnetic 
properties of nanostructures remains to be achieved, de- 
spite the intense effort dedicated to this task. This is 
of intrinsic interest not only for fundamental theoreti- 
cal and experimental physics 1,2 , but also for the practi- 
cal advances that make use of these discoveries to im- 
prove our current technology. In this paper we wish 
to present a contribution to this goal, by outlining our 
finite-temperature theory of the anisotropic spin-spin 
correlations in magnetic metallic monolayers (MLs), de- 
posited on a suitable substrate, and demonstrating its ap- 
plication to Mn MLs on the hexagonal fee (111) surfaces 
of Pd, Pt, Ag and Au. In these systems competing inter- 
actions delicately balance to favour either row-wise anti- 
ferromagnetic (AF) or triangular AF groundstates. Very 
complex unidirectional anisotropies, of Dzyaloshinskii- 
Moriya (DM) type, then conspire to lift the chiral de- 
generacy of the noncollinear state. We extract all this 
information from the analysis of the high-temperature 
paramagnetic susceptibility, which has the spin-spin cor- 
relation function as its main ingredient. 

In itinerant systems there is a mutual feedback between 
the electronic and the magnetic properties. The nature of 
the magnetic state affects the electronic motions, which 
in turn determine the magnetic properties. At T = K, 
first-principles theories such as Density Functional The- 
ory (DFT) provide a quantitative account of the ground- 
state properties 3 . At finite temperature, however, the 
link between electronic structure and magnetism is usu- 
ally broken, by mapping the T = K information to some 
effective magnetic model, which properties are then com- 
puted either by analytical or numerical techniques 4 . 



The theory we present in this paper restores this link, 
by making use of the Disordered Local Moment (DLM) 
picture 5-7 . In the DLM state the relevant degrees of free- 
dom are the local spin quantisation axes associated with 
each lattice site, which can have associated local spin mo- 
ments under favourable conditions. Their orientations 
are assigned according to a specified probability distri- 
bution, and the electrons move through this lattice of 
disordered moments in a mean-field fashion, in keeping 
with the spirit of DFT. The formulation is completely 
general, and can be applied either to bulk or thin films, 
for which the multi-sublattice or multilayer descriptions 
apply 7-10 . 

We focus on the high-temperature paramagnetic (PM) 
state, for which a classical description of the local mo- 
ments is expected to be adequate. Using linear response 
theory, we derive the pair correlation function between lo- 
cal moments, which yields spin-spin correlations similar 
to those obtained in the Random Phase Approximation 5 
(RPA). Previous work along these lines focused on the in- 
terlayer couplings 9 , and we expand on it by also treating 
intralayer couplings, both in real and reciprocal space. 
The fundamental extension of the theory presented here 
is the relativistic formulation, going beyond the isotropic 
correlations by naturally incorporating anisotropic ef- 
fects, which determine the real-space magnetic struc- 
ture and lift degeneracies, such as chirality. In low- 
dimensional systems it is the magnetic anisotropy which 
stabilises the ordered state against spin fluctuations, and 
so it is a crucial ingredient of any theory aimed at de- 
scribing nanostructures. 

The theory is then applied to Mn MLs on the hexag- 
onal (111) surfaces of Pd, Pt, Ag and Au. Our mo- 
tivation stems from recent experiments performed for 
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Mni/Ag(lll), which were interpreted as evidence of a 
triangular AF state 11 , in contrast to the predictions from 
previous electronic structure calculations, which pointed 
to a row-by-row AF state 12,13 , and even ruled out a 
more exotic 3D structure 12 , such as the one predicted 
for Mni/Cu(lll) 14 . By investigating the properties of 
Mn MLs on these four closely related substrates, we un- 
covered a consistent picture of competing magnetic in- 
teractions leading to frustration and different magnetic 
groundstates very close in energy, coupled to very com- 
plex DM-like anisotropics, which favour chiral triangular 
AF states. Our findings are consistent with other pub- 
lished calculations, and suggest that further theoretical 
and experimental work is required to shed light on these 
deceptively simple systems. 

We begin by presenting our new theory of the 
anisotropic spin-spin correlations in magnetic metallic 
MLs. The following section contains the details of our 
numerical calculations, which are presented afterwards. 
We finish by summarising our results and presenting an 
outlook on further work. 



II. THEORETICAL FOUNDATION 



A. Local moments in metals 



Local magnetic moments in metals are an emergent 



property 



5,6,15 



Due to electron-electron interactions, an 



electron in an itinerant system will have an energetic pref- 
erence to move between lattice sites which have an over- 
all spin polarisation more or less aligned with its own 
spin. If the orientations of the local spin polarisation on 
the lattice sites vary more slowly than the timescale of 
the electronic motions, collective degrees of freedom can 
be meaningfully assigned to those lattice sites, which we 
label 'local moments'. The magnitude of these local mo- 
ments will vary in time as fast as the electrons hopping 
to and from each site, but on the timescale of the preces- 
sional motions a well-defined time-averaged magnitude 
is apparent. 

Good local moment systems are those for which the 
magnitude of the magnetic moments is fairly insensitive 
to the local and global magnetic state 16 ' 17 . For these sys- 
tems one can make use of the Rigid Spin Approximation 
(RSA). In DFT language it means that the exchange- 
correlation magnetic field entering the Kohn-Sham equa- 
tions is the same for any orientation of the local moment. 
It is an approximation because there will always be feed- 
back between the electronic and the magnetic properties 
in an itinerant system, however it has been shown in the 
past that for many systems it holds very well, for example 
by performing spin-spiral calculations 12 . We shall make 
use of the RSA in our construction of the DLM state. 
It should be clear, however, that the electronic structure 
in the RSA is still different for different magnetic states, 
and even more so for the DLM state, due to the ensemble 
averages over the local moment orientations. 



B. Statistical mechanics of localised moments 

Let us now develop the meaning of our theory of the 
anisotropic spin-spin correlation function by considering 
the statistical mechanics of a local moment system. The 
Hamiltonian for the system is written as 

JV 

H = U int + 74xt = ^int({e}) - Y, «i ■ ^ (1) 



The relevant degrees of freedom are the orientations of 
the local spin quantisation axes, written as unit vectors 
ii . The first term is the interaction Hamiltonian, describ- 
ing the coupling between the local moments, and is left 
unspecified. The second term is the Zeeman energy due 
to coupling to the external applied fields Hi. 
The free energy is given by 



Z=f[[de k e-f >H = /d{e}< 



-pn 



(2) 



which also defines the partition function, Z. Here 1/(5 — 
kftT, with fee Boltzmann's constant and T the temper- 
ature. Then the local order parameter, m, = (ii), is 
obtained from 



d{e} 



-pn 



d{e}P({e})e, (3) 



dT 



which also introduces the probability distribution, 
P({e}), and the magnitude of the local moments is omit- 
ted. It can be reintroduced in the final expression for the 
susceptibility, if desired. 

The magnetic susceptibility then follows 



drhi 
dH 1 



(4) 



The boldface signifies a 3 x 3 tensor in Cartesian compo- 
nents. This expression defines a spin-spin pair correla- 
tion function, and connects it to the magnetic suscepti- 
bility — it is a special case of the fluctuation-dissipation 
theorem 18 . 

To proceed with the analysis some approximations are 
required. It is already apparent that for the two quanti- 
ties of interest, the local order parameter and the mag- 
netic susceptibility, full knowledge of the complete proba- 
bility distribution is not necessary, and that only reduced 
probability distributions are required: 

(ii) = /de i Pi(ei)e j (5) 

(iiij) = dii dejP 2 (e i ,e j )e l e j (6) 



where, for instance, Pa (e» = Ilfe^ij/defc P({e}). 
Thus a useful approximation may focus on obtaining 
these restricted averages. 
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The other obstacle is the need to handle an arbi- 
trary, possibly highly complicated, interaction Hamilto- 
nian, furnished by the first-principles electronic struc- 
ture. A very successful scheme to address this matter 
is that of variational statistical mechanics, in which a 
trial Hamiltonian, %o({e}), is chosen and its parame- 
ters determined from a variational bound on the free en- 
ergy. This can be obtained from the Feynman-Peierls- 
Bogoliubov inequality 19 , 



[H-Ho), (X) = Jd{e} — X (7) 



so that the ensemble averages are calculated using the 
probability distribution generated by the trial Hamilto- 
nian. The variational parameters are then determined 
by ensuring the equality of the respective restricted 



averages^ 



(n)e z -(n) = (Ho)e,-(n ) (8) 

(WW -(H) = (H )6 t A -{Ho) (9) 



where all orientations are averaged over except those sin- 
gled out, i.e., the orientations on sites i and j are kept 
fixed where stated. 

For simplicity, we choose a site-diagonal trial Hamilto- 
nian, Ho = — J2i hi - ei- Then the so-called Weiss fields 
are given by 



Hi — — / dii ii(n] 



or hi 



ft-^(io) 

am, 



The first definition follows directly from Eq. (8), and is 
computationally tractable. The second one arises from 
the direct minimisation of the trial free energy, and is 
conceptually more useful. 

The local order parameter is then given by the 
Langevin function, 



in. 



coth(/3hi) 



1 



hi = L(fifn) hi, hi = -^ (11) 



By taking its derivative with respect to external applied 
fields {H}, and making use of 



which measures the interactions between the local mo- 
ments. This is the key quantity for our first-principles 
calculations. The statistical mechanics approach we just 
outlined is completely general, and has been applied in 
other order-disorder contexts 20-22 . 

The local susceptibility tensor is given by 



~ S a0 (PM state) (15) 
o 



and can be recognised as the Langevin susceptibility of 
non-interacting local moments. 

Performing a Legendre transform on the free energy, 
replacing the external applied fields by the local order 
parameters as independent variables, we obtain in the 
vicinity of a transition from the PM state into some or- 
dered state 



F{{m^T^T c )*\Y,5™i-Xi} ■ 



8m i 



1,3 

i B J-'8ij. 



SkftT 8u X — \ ■ 8rhj 



(16) 



At sufficiently high temperature this quadratic form is 
positive definite, and so the minimum of the free energy 
is attained for {8m} = {0}, the PM state. The highest 
temperature at which this quadratic form is no longer 
positive definite will then be the transition temperature 
T c into a magnetically ordered state. It is the presence of 
interactions between the local moments, encoded in Sij, 
which drives this transition. 

We proceed by diagonalising this quadratic form in a 
Fourier basis, using 8fhi — il^jdq exp(— i<f ■ Ri) 8m{q) : 



8fhi • ■ 8m j — 



dq 



n 



BZ 



Sm(-q)-S(q)-8m(q) (17) 



where f2ez is the Brillouin zone (BZ) volume/area, solv- 
ing the 3x3 Hermitian eigenvalue problem for each q: 

S(q)-Up(q) = a p (q)u p (q) , p = 1,2,3 (18) 



dfhj drfij dhi 



dHi dhi dHi 



and — ^ = hi 
dhi 



the magnetic susceptibility quickly follows: 

drhi 



dHi 



Xoa k 



Xo,i ' $ik ' Xkj 
k^i 



(12) 



(13) 



where we introduce our spin-spin direct correlation 
function 5 ' 8 , 



;a p = dhf = d 2 (H) 
dmP- dmfdfnE 

J ' 3 



, a,0 = x,y,z (14) 



with eigenvalues cr p {q) and eigenvectors u p (q), which can 
be chosen orthonormal, and define the basis in which 
S (q) is diagonal. Then T c is given by 

3ksT c = maxa p (q) =>■ Sm(q) oc u p {q) (19) 

and the associated eigenvector determines the mag- 
netisation profile associated with the instability. All 
symmetry-related q vectors will have the same eigenval- 
ues, and so in this single-tf picture they are energetically 
degenerate. 

The standard decomposition of a 3 x 3 Cartesian ten- 
sor will be used, and will be found to provide a natural 
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interpretation of the information contained in S^ . It can 
be written as 



S„ = Sh 



si 



Sij + s), 
'ij 



Si 



S ij ~ j [ S V 



l -(s l0 -s%) 



(20) 



The isotropic part is the dominant contribution. From 
<S T (g) = £^cos(<f ■ R\j)S\j it can be seen that it will 
determine the modulation vector q, with corrections from 
the anisotropic terms, if they are strong enough. 

The antisymmetric part of the tensor can be cast into 
the form of an unidirectional, or DM-type, anisotropy: 



• 5m j = Dij ■ (Srhi x Sfhj) 



(21) 



y ! Sfj it is apparent that in 



From S (q) = i V . sin(<f • R 

general this term will give rise to complex eigenvectors, 
and to the notion of chirality. This is demonstrated with 
a simple example in the next section, and with calcula- 
tions for real systems later. 

The remaining uniaxial anisotropy, S^ , determines the 
real space orientation of the local moments, being of 
pseudo-dipolar form. The familiar magnetostatic dipole- 
dipole interaction is of this form, and if it is added to the 
free energy it will be incorporated in this term. 

To put flesh on these statistical mechanics bones 
we will now illustrate this approach for a generalised 
anisotropic Heisenberg model, and compare with more 
familiar treatments. 



C. Spin— spin correlations in the Heisenberg model 

Suppose that the hitherto unspecified interaction 
Hamiltonian is of the form of a generalised anisotropic 
Heisenberg model, 



Hi 



IE 

i,3 



' 3 ij ' &j 



JEE^ 



"/3 £/3 
ij j 



(22) 



hi ct,0 



with a, /3 — x, y, z. Choosing the same trial Hamiltonian, 
we obtain for the Weiss fields 



Hi 



E^, 



(23) 



which reproduces the well-known mean field result. This 
particular choice of trial Hamiltonian will only capture 
the bilinear part of a generalised interaction Hamilto- 
nian, as can be seen from its definition. In this model we 
can identify S^ j7Yp and make immediate use of the 
discussion from the previous section. 

The magnetic susceptibility tensor is then 



Xo,i 



E/ Xo,i 



(24) 



Taking the lattice Fourier transform and inverting to 
solve for the susceptibility, 



X(q)= (XoT'-Jiq) 



PM 



3k B TX-J{q) (25) 



and the paramagnetic limit is shown. Thus the mean field 
theory in an external applied field yields non-local corre- 
lations for the Heisenberg model, which are governed by 
the Heisenberg exchange interactions. 

To finish this section we show an example of what the 
Fourier transformed Cartesian tensor might look like, for 
a particular modulation vector q. Using Eq. (19) to inter- 
pret the tensor, along with u, = Re[exp(— iq ■ Rj) u(q) ], 

/ J|| LD 

J (?) ' u p (q) = (J p {q) u p {q) , J{q)=\-\D J\\ 

V Ji, 

=>Ox{jfi = J|| + D , <7 2 {q) = Jy - D , a 3 (q) = J ± 
-4- Ui(q) = x - iy , ut = cos(g • R4) x - sin(g • Ri) y 
tiziq) = x + iy , u< = cos(g ■ R t ) x + sbx(q ■ R t ) y 
— > uz{q) = z , Ui = cos(g • Ri) z (26) 

This depicts the case of a uniaxial exchange with a DM- 
type contribution equivalent to a DM vector along z. The 
two in-plane eigenvalues have complex eigenvectors, rep- 
resenting two possible chiralities, and their degeneracy is 
removed by the DM anisotropy. These chiral states are 
illustrated in the two bottom panels of Fig. 2. 

We now proceed to derive the magnetic interactions 
from first-principles electronic structure theory. First 
we must construct a description of the high temperature 
paramagnetic state, using the DLM state. Then, by car- 
rying out the same type of analysis, the spin-spin correla- 
tions in the paramagnetic DLM state are derived. These 
take the place of the Heisenberg exchange parameters, 
but can be more general, as will be seen. 



D. The Disordered Local Moment state 

Our theory of the DLM state is based on DFT, in 
the Kohn-Sham picture, which is a single-particle mean- 
field theory. Each electron moves in the self-consistent 
field (SCF) generated by all other electrons. A local spin 
quantisation axis can then be naturally associated with 
each lattice site, according to the average orientation of 
the spin-only magnetic Kohn-Sham potential. The elec- 
tronic properties are obtained from Multiple Scattering 
Theory 23 , and the detailed derivation of the DLM equa- 
tions has been presented before, for the relativistic 7 and 
non-relativistic 5 cases. In what follows only the essential 
definitions and equations will be collected, and we refer 
to the appropriate references for more information. 

The Kohn-Sham potentials for a lattice site i with the 
local spin quantisation axis oriented along define the 
single-site t-matrix, t^ef, e), for given energy e, through 
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the Lippmann-Schwinger equation. Making use of the 
RSA for spherical potentials, the t-matrix for an arbi- 
trary orientation is simply related to that of a given 
reference one, for example along the z-axis, by an appro- 
priate rotation 7 : 



t l {e l ;e) = R{e l )t l {z-e)R\e i ) 



(27) 



Without SOC the spin components factorise, and the t- 
matrix takes on the form 5 



t i (e i ;e)=tf (e) 2 + t t (e) a ■ e l 



(28) 



with a the vector of Pauli spin 1/2 matrices and I the 
identity matrix. It can be shown that then the only pos- 
sible couplings between the orientations of the local mo- 
ments at two different sites is some power of Si • ij , which 
is isotropic in real space. Only by incorporating relativis- 
tic effects can this limitation be overcome 24 , and this is 
the essential new feature of our extension of the theory 
of spin-spin correlations in local moment systems. 

Relativistic effects arise by solving the single-site 
problem using the Dirac instead of the Schrodinger 



equation 



The spin-only magnetic Kohn-Sham po- 



tential, which gives rise to the local spin moment, will 
then also generate a local orbital moment, through the 
SOC effect. The orbital moment is an induced property, 
resulting from the relativistic description, and is not an 
independent order parameter in our theory. 

Propagation between lattice sites is described by mak- 
ing use of a suitably chosen reference system, for which 
the Green's function is known, G^-(e). The scattering 
path operator is then given by 27 

—ij {e)=t i (e) + t_ t (e)J2 (e) Z kj (e) (29) 

and the integrated density of states (DOS) by Lloyd's 
formula 28 , 



N(e) = -ImTr logr(e) 

TT ~ 



(30) 



where the double underlines stand for matrices in site 
and angular momentum labels, and the local moment 
orientations were omitted. 

The resulting Hamiltonian, Ji({e}), for an arbitrary 
set of local moment orientations, {e}, is intractable. 
The theory progresses by invoking a simpler reference 
Hamiltonian, "Ho = — X)i ' an d making use of the 
Feynman-Peierls-Bogoliubov inequality to determine the 
best Weiss fields, as described in the previous sections. 

To carry out the restricted averages the Coherent 
Potential Approximation 29 (CPA) is invoked. In this 
scheme we replace the true disordered medium by an ef- 
fective translationally invariant one, such that the av- 
erage properties of an electron are reproduced. In 
the single-site approximation one requires the effective 
single-site t-matrix, t^e), which is chosen to satisfy 7 



&e i P l (ei)X l {e i ]£) = 



(31) 



and the impurity and excess scattering matrices are 
given, respectively, by 



D i {e i ;s) = 1- [ti (e)-ii (e t ;e) ) f H (e) 



(32) 



and 



X i {e i ;e) = £ L &\ej) 1 (33) 

This leads to the desired expression for the Weiss fields, 
making use of the effective medium quantities: 



-ImTr fdef(e;v) 

7T 



3 

-in 



Jdii ii log = D i (e i ;£) 



(34) 



where f(e; v) is the Fermi-Dirac distribution for energy e 
and chemical potential v. This expression has the form of 
an impurity integrated DOS, and in fact is derived from 
Lloyd's formula 7 , Eq. (30). 

Now that we presented our construction of the DLM, 
we can proceed and perform a linear response analysis 
of how the DLM state responds to an inhomogeneous in- 
finitesimal external magnetic field applied in each site. 
This procedure will generate our desired spin-spin cor- 
relation functions, which contain information about the 
coupling between magnetic moments, determined by the 
self-consistent effective medium. 



E. Linear response and magnetic susceptibility 

We begin by formally introducing an infinitesimal ex- 
ternal magnetic field on each site, 5 Hi. The Weiss field 
on site i, hi, will then be modified as hi — > hi +SHi+Shi. 
The small correction Shi is due to the adjustment of the 
self-consistent effective medium in the presence of the 
infinitesimal external fields. In the absence of this feed- 
back the magnetic response would be that of a system of 
non-interacting Langevin spins, as we will now show. 

The response of the local magnetisation to the external 
field is by definition the magnetic susceptibility. In our 
linear response approach it is given by 



Sfhi 



Xa,t 



s i:i x 



dhj 



X Q ,-W tl (35) 



with the magnetic susceptibility of non-interacting mo- 
ments defined as before. All probability distributions 
and effective medium quantities are those of the reference 
state, due to the linear response approach being used. 

The coupling between the local moments is then con- 
tained in the response of the local Weiss fields to the 
infinitesimal external fields. Using the definition we ob- 
tain, after some algebra, 



Tr 



47T 



su 



■/ k 



sm 



-hi 



(36) 



G 



where the energy integration, taking the imaginary part 
and dividing by ir is being implied along with the trace. 

A second equation can be obtained from the equation 
that defines the DLM effective medium, Eq. (31), by ob- 
taining the first-order correction to the effective medium 
i-matrix. Introducing the expansion 



5Ht 



(37) 



which defines a 'full' vertex, Ay, the following equation 
of motion can be written: 



Ay — A l 5%j + 



i) JZ^ife Afcj Tki X-i(Zi) 

(38) 



by defining the 'bare' vertex 



Aj = dii Pi(ii) (e* - rhi) 2Ci(e») 



(39) 



This gives for the magnetic susceptibility the expected 
expression 



Xij — Xo,i + g 2.*/ ik ' Xkj 

k 



(40) 



where the spin-spin correlation function is given by 
Stf = -ImTr Jdef(e;v)A-(e)Y^t tk (e)f kj (e)f kl (e) 

(41) 



k=£i 



with A t (e) = — I d£i £i Xiief, 



If relativistic effects are not included in the theory, then 
this expression becomes independent of the Cartesian in- 
dices, and thus isotropic in real space. This is a con- 
sequence of the previous discussion on the form of the 
single site i-matrices with and without SOC effects. 

From the corresponding expressions, Eqs. (38) and 
(39), it can be seen that the 'bare' vertex is the direct 
response to the external applied fields, while the 'full' 
vertex also incorporates fluctuation terms, defined by the 
angular averages in Eq. (38). For the systems being 
studied in this paper these are small corrections, and so 
we can approximate the 'full' vertex by its first iteration, 
by replacing it by the 'bare' vertex on the right-hand- 
side of Eq. (38). 

Focusing on the high temperature PM state, and make 
use of its high symmetry to take the lattice Fourier trans- 
form of Eq. (40) , and inverting to solve for the suscepti- 
bility: 



X (q)=V 2 3k B TI-S(q) 



(42) 



where the magnitude of the local moments in the PM 
state, //, was reintroduced. This shows that the spin- 
spin correlation function plays a similar role to that of 



the Fourier transform of the exchange constants in the fa- 
miliar Heisenberg model. The analysis that was carried 
out for that case can then be mapped to the present de- 
scription, and be used to interpret our results. It should 
be noted, however, that the spin-spin correlation func- 
tion is derived from the electronic structure in the DLM 
state, and incorporates fluctuation effects. 
The core of Eq. (41) is the convolution 



T-ik Afcj Tki 



dk 



BZ 



f(k+q)\(q)m (43) 



which shows that Fermi surface nesting effects are ac- 
counted for, through the product of the scattering path 
operators. Thus RKKY behaviour may occur 30 . 

We will now present the results of our calculations for 
Mni/X(lll), X = Pd, Pt, Ag and Au. This choice of 
closely related substrates will highlight the role of hy- 
bridisation with the substrate in determining the pre- 
dicted magnetic instability, as well as the role of varying 
SOC on the anisotropic effects. 



III. COMPUTATIONAL DETAILS 

Our calculations are carried out in the framework 
of the Screened Korringa-Kohn-Rostoker 31 (SKKR) 
method for layered systems. The relativistic formalism 
was adopted, in the Local Spin Density Approximation 
(LSDA), using the atomic sphere approximation. The 
energy integrations were performed on an asymmetric 
semicircular contour in the upper half of the complex 
plane using 24 points. 

The self-consistent paramagnetic DLM potentials were 
generated in the scalar-relativistic approximation, using 
200 k-points in the irreducible wedge of the 2D Brillouin 
zone, which has C§ v symmetry. These potentials were 
then used in the fully relativistic linear response calcu- 
lations, with 3721 k-points in the full 2D Brillouin zone 
ensuring all results for the spin-spin correlation function 
are converged within 0.1 meV. This is possible as only 
complex energies with finite imaginary part are used in 
the energy integration. 

The geometry used in our calculations was as follows. 
The interface region was composed of eight layers of the 
transition metal substrate, the Mn ML, and three layers 
of empty spheres. This was then matched to the semi- 
infinite bulk substrate and to the semi-infinite vacuum 
region, so that no slab approximation was used. 

No attempt was made at determining the equilibrium 
geometries for the four systems studied. However, the in- 
fluence of different interlayer spacings between the mag- 
netic monolayer and the substrate was investigated for 
the cases of Mn/Ag and Mn/Au, and it was found that 
there was no qualitative change in the magnetic proper- 
ties, up to 15% inward relaxation. The dependence of 
the magnetic properties on the in-plane lattice constant 
was addressed in a previous study 32 , and for the range 
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FIG. 1. The isotropic part of the S(q) tensor along high 
symmetry lines of the 2D BZ. The maxima signal the favoured 
magnetic instability, which correspond to the row-wise AF 
(Ag, Au) or triangular Neel states (Pd, Pt), see also Fig. 2. 



relevant to our systems the same conclusion can be made. 
On the other hand, the relaxed geometrical parameters 
depend on the magnetic state and other approximations, 
such as the exchange-correlation functional and the han- 
dling of the charge density, so a definite statement is 
precluded. We thus progress with a simple model of the 
geometric structure, which we now describe. 

To aid comparison, the average of the experimental 
lattice constants for Pd and Pt was used for both Pd and 
Pt (a = 3.905 A), and likewise for Ag and Au (a = 4.085 
A) . This means that the in-plane lattice constant for the 
Mn ML is the same for the Pd and Pt substrates, and also 
for the Ag and Au substrates. To estimate the inward 
relaxation, it was assumed that Mn grows on Cu(lll) 
without any significant relaxation, thus defining the Mn 
ML hard sphere radius. By considering how these hard 
spheres stack on top of the larger ones corresponding 
to the substrate, the following relation was obtained for 
the ratio between the esti mated and the ideal dist ance 
between planes: d/d = \/3/8 (1 + ac u /ax) 2 — 1/2. We 
then round the results and use 5% inward relaxation for 
the Pd and Pt substrates, and 10% for Ag and Au. 




FIG. 2. Magnetic structures associated with high-symmetry 
points on the 2D BZ. Top: row-wise AF state associated with 
the M point. Middle and bottom: triangular Neel states of 
opposite chiralities — see discussion in text. Spheres show the 
positions of the magnetic atoms. Crosses mark the positions 
of the nearest substrate atoms. 'Up' and 'down' triangles are 
shown, with respect to the nearest substrate atoms. 



A. Instabilities of the paramagnetic state 

First we discuss the isotropic part of the tensor, «S I ((?). 
This can be conveniently plotted along high-symmetry 
lines in the 2D Brillouin zone, and can be used as a guide 
to the favoured magnetic instability, given by the modu- 
lation vector associated with its maximum value, see dis- 
cussion following Eq. (16). For the Pd and Pt substrates, 
the favoured magnetic instability is associated with the 
K and K' points, which corresponds to triangular Neel 
states, while for the Ag substrate the favoured instabil- 
ity is associated with the M point, which corresponds to 
the row-wise AF state. The case with the Au substrate 
appears to show an approximate degeneracy between the 



IV. ANISOTROPIC SPIN-SPIN 
CORRELATIONS IN Mn MONOLAYERS 

We begin by presenting our results for the anisotropic 
spin-spin correlations in Mn 1 /X(lll), with X = Pd, Pt, 
Ag and Au, and then proceed to a discussion and com- 
parison with other theoretical and experimental works. 
The same decomposition that was used in Eq. (20) will 
be used to extract the information from the spin-spin 
correlation function. 



TABLE I. Calculated real-space isotropic interactions, Slj, 
and magnitude of DM vectors, |.Dy|, for first and second NNs. 







Mn/Pd 


Mn/Pt 


Mn/Ag 


Mn/Au 




First shell 
Second shell 


-53.33 
-2.50 


-55.49 
-2.72 


-18.51 
-5.91 


-22.92 
-3.22 


\Dij\ 


First shell 


1.19 


4.41 


0.29 


0.34 


Second shell 


0.34 


0.27 


0.00 


0.28 



K and M point-type of AF. 

The high symmetry magnetic states are illustrated in 
Fig. 2. Considering only the information in S l (q) , as K 
is mapped into K' by a mirror symmetry, they are equiv- 
alent, and degenerate in energy. Nothing can be said 
about the real space orientation of the local moments, 
only about the relative angles between them, as the ten- 
sor is isotropic in real space. 

Looking at the real-space values of the tensor, see Ta- 
ble I (a negative sign means an AF coupling), it can be 
seen that the contribution from the nearest-neighbours 
(NNs) is dominant for Pd and Pt and, though smaller in 
magnitude, also for Ag and Au. These interactions leads 
to the competition between magnetic states shown in Fig. 
1, for the Au substrate. This competition effect is similar 
to the behaviour of the J1/J2 isotropic Heisenberg model 
in a triangular lattice 12 . 



B. Real space magnetic structure 

The anisotropic part of S (q) is small compared to the 
isotropic part, but it sets the real space orientation of the 
spins. The largest eigenvalue of S(q) has local extrema 
at the same positions as <S (q), see Fig 1, so the magnetic 
instabilities remain unchanged. For the systems studied, 
the symmetric anisotropy S s (q) favours an in-plane spin 
configuration. 

The antisymmetric part of the tensor, S A (q), can be 
cast in the form of a DM vector, see Eq. (21). The size 
of this DM vector is the same for all atoms in the same 
lattice shell, and is plotted as a function of the distance 
between the atoms in Fig. 3. As expected, it is largest 
for the first shell, and also greatly enhanced in the Pd 
and Pt substrates. The case of the Au substrate is most 
interesting, as there are several lattice shells with com- 
parable DM vectors (see Fig. 4), leading to a complex 
unidirectional anisotropy field. In the S(q) picture, this 
anisotropy leads to a splitting of the eigenvalues, as ex- 
emplified in Eq. (26). 

If the triangular Neel state is the favoured magnetic 
instability, the presence of the DM anisotropy will select 
one of the two possible chiralities. This is the case for 
the Pd and Pt substrates. We define a chirality vector in 
the following way 33 : 



£ = Smi x 8rri2 + 5rfi2 X 5rn 3 + (Sm 3 x 5rhi 



(44) 



Here the numbers label the magnetic atoms going an- 
ticlockwise around each triangle of NNs in the lattice. 
There are two types of triangles in the magnetic layer: 
those that do encircle an atom of the nearest substrate 
layer, which will be called 'up', and those which do not, 
called 'down'. These two types of triangles are marked in 
Fig. 2. For collinear magnetic states the chirality vector 
is obviously zero (see top panel of Fig. 2). 

Referring to the middle panel of Fig. 2, for all 'up' tri- 
angles £ oc — z (the 'spins' rotate in a clockwise fashion), 




2 3 4 5 6 

Distance (units of 2D lattice constant) 



FIG. 3. The absolute value of the DM vector against dis- 
tance. In the Pd and Pt substrates the DM vector for the 
first shell of neighbouring atoms is greatly enhanced, see Ta- 
ble I. The existence of sizable DM vectors over several shells 
is remarkable. 



-A. — V— f — I 1 — f 



FIG. 4. The DM anisotropy field for Mn/Au. Vectors of the 
same colour have the same magnitude, and belong to the same 
lattice shell, representing the interaction with the central site. 
Note that z-component is small for most DM vectors, and so 
is not very apparent in this 3D plot. 



and all 'down' ones £ oc z (the 'spins' rotate in an anti- 
clockwise fashion) . This pattern is reversed for the mag- 
netic structure depicted in the bottom panel. The net 
chirality is thus zero for both magnetic states. However, 
these two chirality patterns are clearly distinguishable, 
as the 'up' and 'down' triangles are inequivalent. This is 
encoded in the C% v symmetry of the system, and allows 
the lifting of the chiral degeneracy. 

It can be shown straightforwardly that the high sym- 
metry of the triangular Neel state dictates that the en- 
ergy difference between the two chiral states is set by 
the z -component of the DM interaction 34 . For the sys- 
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tems studied, the DM vectors are mostly in-plane, thus 
their z-component is small, as can be seen for example 
in Fig. 4. This means the lifting of the chiral degeneracy 
is smaller than would be anticipated, from considerations 
based on the length of the DM vectors. The energy dif- 
ference between the two magnetic pattenrs depicted in 
the lower panels of Fig. 2 is of the order of 1 meV in 
all cases studied. It is not large enough to favour the 
triangular Neel state, when the isotropic part of the ten- 
sor favours the row-wise AF state, as in the Mn/Ag and 
Mn/Au cases. If the dominant isotropic interactions were 
to favour another noncollinear state with less symmetry, 
as was demonstrated for Mni/W(001), then this splitting 
would be more significant 35 . 



Ag in Mn/Ag 
Mn in Mn/Ag 
Au in Mn/Au 
Mn in Mn/Au 




-4 -2 

Energy from the Fermi level feV) 



C. Connecting magnetism and electronic structure 

We now proceed to establish a link between the elec- 
tronic structure and the just described magnetic proper- 
ties. 

The influence of the substrate on the electronic struc- 
ture of the Mn ML can be readily understood from the 
layer-resolved DOS, in the DLM state, shown in Fig. 
5. The Mn minority states are just slightly affected by 
changing the substrate, as can be seen be comparing all 
four cases. The majority states, however, show a marked 
difference. For the Ag and Au substrates these states are 
fairly narrow, while for the Pd and Pt substrates they 
are much wider in energy, and show strong signs of hy- 
bridisation. This can be explained by the position of the 
substrate d-states. Ag and Au have filled d-bands, lower 
in energy than the Mn majority states, and so there is 
little hybridisation. On the other hand, Pd and Pt have 
partially filled d-bands, which extend up to the Fermi 
energy, so there is much stronger hybridisation with the 
Mn majority states. The strength of the hybridisation 
with the substrate alo affects the size of the DLM spin 
moments: for Mn/Ag and Mn/Au it is about 3.9 /ie, in 
agreement with previous calculations, while for Mn/Pd 
and Mn/Pt it drops to about 3.6 [1b- 

A simple picture can be invoked to synthesise these 
arguments: on the Ag and Au substrates, the magnetic 
properties of the Mn ML are quite similar to those of 
a conceptual free standing ML, as found before in other 
electronic structure calculations. This explains the sim- 
ilarity between the S\, ; tensors for the Ag and Au sub- 
strates, see Table I. The Pd and Pt substrates then rep- 
resent a strong enhancement of the NN interactions due 
to hybridisation with the substrate. 

Referring also to Fig. 3, it can be seen that the 
anisotropic effects scale with the atomic number of the 
substrate, due to the SOC effect, as expected. The DM 
vectors on the first shell, however, show a marked differ- 
ence: they are much larger in the Pd and Pt substrates 
than in the Ag and Au substrates, regardless of the sim- 
ilarity between atomic numbers. This reveals another 
ingredient to a large anisotropic effect: strong hybridis- 
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FIG. 5. The layer-resolved DOS for Mnl/X(lll), with X 
= Pd, Pt, Ag and Au. Only the DOS for the Mn and the 
nearest substrate layer is shown. 



ation with the substrate. This is why the first shell DM 
vector in Mn/Pd is three times larger than the one in 
Mn/Au, despite the differences in atomic number, and 
so in SOC strength. 



D. Mni/Ag(lll): comparison and discussion of 
previous theoretical and experimental work 

Now that we presented our results for the four systems 
studied, we wish to focus on a more detailed comparison 
for the most studied case, that of Mni/Ag(lll). 

It was found experimentally that the magnetic ground 
state was the triangular Neel state 11 , in contrast with the 
row-wise AF state predicted by our theory, and also by 
previous electronic structure calculations 12,13 . We sug- 
gest two explanations for the disagreement between the- 
ory and experiment. Both are based on the delicate bal- 
ance found between the two types of AF: 

• There may be some buried Mn atoms underneath the 
studied Mn triangular islands. This would enhance 
the hybridisation to the substrate, which has a very 
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strong influence on the strength of the NN couplings, 
as found for the Pd and Pt substrates, and thus favours 
the triangular AF state. 

• There may be an inadequate treatment of the 
exchange-correlation effects in the Mn ML. Although 
this might be a small correction, the Mni/Ag(lll) 
system, as described by the LSDA, has a small en- 
ergy difference between the two types of AF states just 
discussed. The electronic structure changes that may 
be brought about by an improved treatment might be 
enough to tilt the preference from the row-wise to the 
triangular AF state. 

The experimental data on Mni/Ag(lll) also suggests 
that similarly oriented magnetic islands with respect to 
the substrate have different magnetic domains, and so 
different magnetic anisotropies. We propose that this 
might be due to two different chirality patterns of the 
triangular AF state. As explained by the experimental- 
ists, similarly oriented islands may grow on the substrate 
in an fee or hep-like stacking. Going back to our dis- 
cussion on chirality, the difference between the two cases 
is whether the nearest substrate atoms are in the cen- 
tre of the 'up' or 'down' triangles. Assuming that the 
'spins' around the substrate atoms always have the same 
chirality vector, the two different ways in which the Mn 
atoms can stack on the substrate lead naturally to two 
distinct chirality patterns, and so two different magnetic 
domains. 



V. CONCLUSIONS 

We presented a new theory of the anisotropic spin-spin 
correlations in good local moment systems, and demon- 
strated its application to Mni/X(lll), with X = Pd, Pt, 
Ag and Au. The kind of magnetic structure for which 
the PM state of each system is unstable was identified. 
Mn/Pd and Mn/Pt should order in the triangular Neel 
state, while Mn/Ag and Mn/Au order in the row-wise 
AF state. For the last two systems it was found that the 
energy difference between the two types of AF is quite 



small. The symmetric part of the anisotropy dictates 
that these should be in-plane magnetic configurations, 
while the antisymmetric, or DM-type, anisotropy lifts 
the chiral degeneracy of the triangular Neel state. 

For Mn 1 /Ag(lll), which is the most studied system 
in the literature, we found agreement with the findings 
of previous theoretical work, while experimental results 
point to a different magnetic state, namely the triangular 
Neel state. We suggested two resolutions for this incon- 
sistency. Moreover, we proposed that the two types of 
magnetic domains found experimentally are manifesta- 
tions of the chirality of the triangular AF state. 

We also wish to highlight our observation that the 
anisotropic effects induced by the substrate do not sim- 
ply scale with its atomic number. The strength of the 
hybridisation with the substrate may play a dominant 
role, and so a careful choice of substrate might yield sim- 
ilar anisotropic properties using lighter elements, thus 
dispensing the use of their heavier counterparts, which 
may be rarer and more expensive. 

The theory presented here is more general than the ap- 
plications chosen to illustrate it. In the interest of sim- 
plicity the formalism was presented for a single magnetic 
lattice or layer, but can be straightforwardly extended 
to multisublattice or multilayered systems. Our investi- 
gations are now progressing in this direction, with the 
Co/IrMn exchange-bias system, a spin valve element, as 
its driving application. 
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